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Abstract 



We use perturbative QCD, beyond the leading InQ^ approximation, to show how 
measurements of diffractive J/il) production at HERA can provide a sensitive probe of 
the gluon density of the proton at small values of Bjorken x. We estimate both the effect of 
the relativistic motion of the c and c within the J / ip and of the rescattering of the cc quark 
pair on the proton. We find that the available data for diffractive J/ip photoproduction 
can discriminate between the gluon distributions of the most recent sets of partons. 



1. Introduction 

The observation of high energy diffractive J / ijj photo- or electroproduction, 7*-*-'^ J/'ip P, 
offers a unique opportunity to measure the gluon density in the proton at low x. Indeed, for 
sufficiently high 7p centre-of-mass energy W, perturbative QCD can be used to express the 
cross section for this, essentially elastic[|, process in terms of the square of the gluon density. 
To lowest order the 7*p J/ip p amplitude can be factored into the product of the 7 — > cc 
transition, the scattering of the cc system on the proton via (colourless) two-gluon exchange, 
and finally the formation of the J/ ip from the outgoing cc pair. The crucial observation is that 
at high W the scattering on the proton occurs over a much shorter timescale than the 7 ^ cc 
fluctuation or the J/ip formation times, see Fig. 1. 
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Figure 1: Schematic diagram for high energy diffractive J /if) photoproduction. The factorized 
form follows since, in the proton rest frame, the formation time Tf ~ 2E^/{Q'^ + M|) is much 
greater than the interaction time Tint — R where R is the radius of the proton. 

Moreover this two-gluon exchange amplitude can be shown to be directly proportional to the 
2 

gluon density xg{x, Q ) with 

= (g2 + Mj)/4, X = AQ'/W'. (1) 

Q'^ is the virtuality of the photon and is the rest mass of the J/tp. To be explicit, the 
lowest-order formula is|[l[] 

The derivation of the result is sketched in section 2(a). An analogous formula to (^ was 
presented by Brodsky et al. but only for longitudinally polarized vector mesonsQ. An earlier 

^Elastic in the sense that the photon and J/^p have the same quantum numbers. 

^It has been pointed out j|] that a factor of 4 should be included in the numerator of the formula in ref. 
Then the results of refs. and |^ agree. 
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work PI presented a "soft" Pomeron treatment of diffractive vector meson production, but 
in this case the connection to the gluon distribution g{x, Q ) was not made. Here we note 
that the heavy meson mass should ensure that perturbative QCD can be apphed even 
in the photoproduction hmit, = 0. The last term in (0) allows for electroproduction via 
longitudinally polarised virtual photons with ai/cTT ~ Q"^ /M"^. The result @} assumes a non- 
relativistic wave function for the J/ifj with the c and c having momenta . Eq. (H) is derived 
assuming the leading InQ^ approximation in the integral d'^k over the gluon loop in Fig. 1. 
When, in section 2, we improve on this approximation we find from the explicit form of the 
integral that typically ^ Q . It is this which specifies the scale of as in (H). 

We may use the available measurements of the 7p J / ipp production cross section (together 
with the observed J/ip diffractive slope b = 4.5 GeV^^) to give a first estimate the of gluon 
density. The results are shown in Fig. 2. 



Figure 2: The points show the lowest-order estimate of the gluon density obtained from the 
available high energy ai^jp — > J/'ip p) data [5,6] using eq. (2). For comparison we also show the 
gluon of the GRV parton set [7], and of the MRS(A', G) partons of ref.[8] evolved back to 
= 2.5 GeV2. 
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This figure is simply to illustrate the typical precision and the kinematic range of the gluon 
density that is probed by the J/ip data. We emphasize that it is a first look at the gluon and 
that the errors reflect only those of experiment. The curves correspond to gluons obtained 
from the latest sets of parton distributions. The spread of the predictions demonstrates the 
potential value of this process as a measure of the gluon. We will show that a comparison of 
the shape of the curves with the data is more reliable than that of the normalisation. Thus, at 
this preliminary stage, we see that the shape of the J /if) data favour the MRS (A') gluon. The 
purpose of this paper is to refine the theory so that a meaningful probe of the gluon distribution 
can be obtained. However, it is already clear that diffractive J / ip production at HERA will offer 
a unique, precise probe of the gluon, xg {x,Q ), in the critical low x, Q > 2.5 GeV^ region, 
provided that we can improve on the validity of (||) . 

In section 2 we scrutinise the approximations used to derive and, more important, we 
implement corrections so that quantitative information on the gluon density can be obtained. 
Then, in section 3, we illustrate the discriminatory power of the diffractive J/tp production 
data by comparing with the predictions of the cross section obtained from the gluon densities 
of the latest parton sets. 

The gluon density at low x has so far been constrained by the slope of the deep inelastic 
structure function dF2/dlnQ'^ . We find that J/ip photoproduction is a much more sensitive 
measure. This is only to be expected since the J/ip photoproduction cross section depends 
quadratically on the gluon {a ~ g"^), whereas in deep inelastic scattering we have a linear 
dependence and then only on the derivative of a cross section {da ~ g). 

2. Improved formula for diffractive J/i/j production 
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Figure 3: Lowest-order perturbative QCD diagrams for diffractive J/ip production. The 
particle four-momenta are shown; z and 1—z denote the fractions of the longitudinal momentum 
of the photon that are carried by the c, c quarks. The small crosses indicate which quarks may 
be regarded as essentially on-mass-shell. 



The amplitude for 7'^*^^ — > J/tp p is obtained, in lowest-order perturbative QCD, from 
the sum of the two diagrams shown in Fig. 3. Formula (0) is based on the leading InQ^ 
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approximation in which we assume that the gluon transverse momenta kx satisfies ^ Q . 
Moreover the non-relativistic form is taken for the J/ip wave function 

ij'{z,>,T) = 6^'\kt)6{z - i) (3) 

where z = ^ + K^^/qJ. In other words, the c and c quarks are not allowed to have any Fermi 
momentum inside the J/ip, i.e. k is set to zero. We discuss these, and other, approximations in 
turn below. In particular we compute the corrections which should be applied before confronting 
(|^) with the data. 

(a) Beyond the leading InQ^ approximation; inclusion of the gluon kx 



Since we are primarily concerned with diffractive J/ip photoproduction at small x we work 
in the (leading order) In 1/x approximation, and retain the full dependence and not just the 
leading InQ^ component. We must therefore express the cross section in terms of an integral 
over the (square of the) unintegrated gluon distribution f{x,k^), and so retain the explicit 
gluon kr dependence. 

We first evaluate the gluon loops in the Feynman diagrams shown in Fig. 3. It is convenient 
to perform the loop integrations in terms of Sudakov variables. That is the various particle 
four momenta are decomposed in the form 

ki = aiq' + (3ip' + kiT (4) 

where p' and q' are the proton and photon light-like momenta: p'^ = g'^ = and W"^ = 2p'.q'. 
In particular 

p = p + apq', q = q' + (3^p' (5) 

with ap = rrip/W^ and (3^ = —Q^/W"^. Within the non-relativistic approximation the quarks 
of momenta ± /t are almost on- mass-shell. 

The integration over the gluon longitudinal momentum puts, in the first diagram, the upper 
quark with momentum h = ^q'^ + K — k' on-shell, leaving only the quark propagator (r^ — m^)~^ 
to be integrated over in the gluon kr integration. To express the propagator in terms of k^ we 
first note that 

r'^ = [q- hf = q^ - 2q.h + m^. (6) 
Using the Sudakov decomposition 

h = ahq' + PhP' - k'rp (7) 

we obtain 

2q.h = {f3^ah + f3h)W^ = -zQ^ + {ml + k'S)/z, (8) 

where we have neglected kt- Moreover in the non-relativistic approximation z = and so from 
d) and (D we find 

- ml = -2g^ - 24 (9) 
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where Q = (Q^ + M|)/4, with = —q^, as usual. Here we take ~ 4m^. In Fig. 3(b) 
only the quark with momentum r' is off-shell. In analogy to the derivation of (^ we find 

r'2 - ml = -2Q\ 

Thus the forward scattering amplitude for diffractive J/ip production from transversely 
polarized photon is 

r dk"^ / 1 1 \ 

^^.4.U/,../^(^-^^-^)GWe.«. (10) 

where colour factors give rise to the opposite sign of the two diagrams, and where the amplitude 
is defined by 

The constant gj specifies the cc coupling to the J/ip and Cc is the charge of the c quark. The 
coupling gj may be determined from the width Fee of the J/ip ^ e^e~ decay, a process which 
is described by the same cc quark loop structure. We have 

2 2 Fge-^W^ 



Via 



where = | Vna. The function G{}z) specifies the probability of finding the gluons in the 
proton. In the simplest three valence quark model 

G = ^ ^ 3 (11) 

where | is the colour factor. 
In the realistic case 

G{k) = /bfkl(x,4) (12) 

where /bfkl is the gluon density unintegrated over kx that satisfies the BFKL equation, an 
equation which effectively resums the leading as ln(l/x) contributions. To gain insight into this 
identification, and to specify the value of x, it is sufficient to study the one-rung contribution 
sketched in Fig. 4. 

To determine the momentum fractions x and x' of the gluons in Fig. 4 we note that 

X + X = (3^ - f]^ = — ^^^^ (13) 

X' = /?out - An = (14) 

where the A are the fractions of the p' momentum carried by the various particles i (see (P) 
and where for simplicity we neglect in the quark loop. 
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Figure 4: One of the one-rung diagrams which give the a;s'ln(l/x) contribution to the BFKL 
gluon, /BFKL(a;, k"^)- The emitted gluons have proton momentum fractions x and x' with \x'\ -C 

X. 

Thus, we see that -C x. The entire^ origin of the x dependence of Fig. 4 comes from the 
integral over the rapidity of the s-channel gluon 

' ^ (15) 

Xg 

where the limitQ comes from = xi — x ~ xi > x. This is the origin of the BFKL ln(l/x) 
contribution^. Diagrams with further gluon rungs give agln"(l/x) contributions which are all 
resummed by the BFKL equation for unintegrated gluon density /BFKL(a;, k^), with 

X ^ {Q' + M^)/W^ = 4Q^/W^. 

To relate /bfkl to the conventional gluon density, which satisfies GLAP evolution, we must 
integrate over k^ 



xg{x,Q'^) = I ^ /bfkl (a;, /ct), (16) 



rQ' dkl 

xg{x,Q^} = J 
and so the 7'^*''^ J/ i^P forward amplitude (|T0|) becomes 



■ [°° dkr d{xg{x,kV)) 

^ ^ Jo 2Q\q' + 4) dkl ^ ' 

The integral converges for k^ ^ Q , but the gluon distribution X5'(x, k'^) is not known as 
^ 0. We therefore write 



A ~ i27r^ CcQj as 



xg{x, Ql) f°° dk^ dxg{x, 



+ 



JQI Q'iO' + kl) dkl 



(18) 



■^Recall that spin-one t channel gluon exchange leads to an energy independent (that is x independent) cross 
section. 

^In general the limit is Xg > max{x, x'}, but in our case x 3> \x'\. 

^The virtual diagrams which lead to the Reggeisation of the t channel gluons also have the same ln(l/a;) 
structure. 



where in the first term in the brackets we have neglected the effects of in comparison with 
Q . Eq. (|1^) allows the effects of the gluon kr to be estimated. Eq. (18) is true to 0{Q1/Q ); we 
investigate below the sensitivity of our results to variation of the choice of Qg. The modification 
of (^) becomes apparent when we use ( |T^) to calculate the forward production cross section. 
We have 

(19) 



Tee Ml vr^ , 
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where the [. . .] contain the entry shown in brackets in ([18|). In fact the accuracy of formula 
(19) is even better than the BFKL approximation since the main part of the corrections may 
be hidden inside the experimentally determined values of g{x, Q"^) and Tee- To be more precise, 
if we neglect the tmin effects of section 2(b), eq. (17) is valid for a gluon distribution obtained 
from any evolution equation. 

(b) Discussion of t^^^ effects 

We have tacitly assumed that we are considering a forward "elastic" scattering amplitude 
with t = 0. However, for 7p J/ipp the minimum value of \t\ is 

2_2 




t^in = — TTT^ m, ~ x'mi. (20) 



This result is evident from (|T^); we have to transfer longitudinal momentum through the 
t-channel two-gluon exchange. We have already checked that the difference between the mo- 
mentum fractions x and x' of the two gluons does not affect the BFKL identification ([T^) of 
G; it contributes beyond the leading ln(l/a;) approximation. For GLAP we could, in principle, 
recalculate Fig. 4 with Altarelli-Parisi kernels [^] which take into account the difference x ^ x' . 
Alternatively we can estimate the effect from the analytic structure of A in the complex t plane, 
where the nearest singularity is the I-k threshold at to = 4m^. Extrapolation from tmin to t = 0, 
where the identification xg{x^Q ) is true, gives at most a correction of order tmin/i^o- In feet 
the 2tt singularity is weak, and to ~ is a more representative value |iy]. In any case we 
expect the tmin effects to be small for a; < 0.1. 

(c) Relativistic effects in the J/ip wave function 

Here we discuss the effect of the motion of the quarks within the J/ip. This Fermi motion has 



been considered in ref. [O] where it was concluded that it will lead to a significant suppression 
of the cross section for diffractive production. Only longitudinal vector meson production was 
studied there, whereas for photoproduction we have to deal with transversely polarised J/ip 
mesons. 



7 



To estimate these relativistic effects we combine knowledge of the J/ip wave function 
ip-^lKx, z) with the structure of the couphng of the photon to the cc pair. For transverse 
and longitudinally polarised photons the amplitudes are of the form 



J/ipp) oc — 



(Q2 + 4)J 



xg{x,Q'^ + K^) '4>-^{kt-,z) 
X (Put dz 



(21) 



2z{i - z) y/g^ (g^ - 4) 
(g2 + 4)3 



xg{x,Q'^ + Kt)^1)\kt,z) d^KT dz (22) 



where = z{l — z)Q'^ +m1 and x = {Q'^ + M^)/W'^. These expressions may be compared with 
the simplified non-relativistic approximation in which we take z = ^ and kt = 0, see (^. Thus, 
if we have knowledge of the J/ ip wave function we can compute the correction factor F"^ arising 
from the use of a more realistic J/tp wave function. In particular, for J/ip photopro duct ion we 



have from (^1|) 



/ (1 + 4) ^ xg{x, "^0(1 + ■'^r)) ^■^('^T, z) dvj 



with 4 



2^2 

3*^ 



r ,J( ^. 2 (23) 

xg{x,m^) j f-^{KT, z) ''-^ 

K^/ml where v is the velocity of the quarks in the J /ip rest frame. In writing 
(23) we have set 2; = | in the numerator, consistent with the approximation of neglecting 
contributions of (9(4)- 

If we take a Gaussian form for the J /if) wave function, tl)"^ = 74exp(— a 4/^0)5 then 
f^) = l/2a. Estimates of (f^) obtained from studies of charmonium |13| vary in the 



(4) 



range 0.12 < (v^) < 0.25, with for example a recent lattice calculation giving (f ^) values in 
the interval 0.18 to 0.12 as rric increases from 1.45 to 1.85 GeV. If we evaluate (23) using these 
ranges of (f ^) and rric for each of the three recent sets of partons [GRV, MRS(A', G)] then we 
find that the suppression F"^ of the cross section for diffractive J/ ijj photoproduction at HERA 
lies in the interval 

0.4 < F^ < 0.6. (24) 

The suppression due to the factor (1 + 4)^^ is partly compensated by the larger scale at 
which the gluon is evaluated. An independent study of the relativistic corrections to J/ip 
photoproduction has been made by Jung et al. ITH]. Using the gluonic and leptonic widths of 



^Here we neglect kr in comparison with kt- It is straightforward to derive formulae in the presence of both 
kx and kt- To simplify the presentation we choose to investigate these corrections in turn and hence impose 
the leading logQ^ approximation in eqs. (21) and (22) [1,2]. Correlated effects will occur, but at the level of 
other uncertainties. 

^Eqs. (21) and (22) are derived assuming the form of the J/4> light cone wave function given in (2.22) of ref. 
[2] . The derivation is subtle and will be presented elsewhere in a more general study of diffractive vector meson 
production[|l2t. 
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the J /ip they estimate {v"^) = 0.16 with = 1.43 GeV, values which are quite compatible with 
the above estimates. The correction factor F"^ evaluated using these values of (f ^) and rric is 
shown in Table 1 over the range of x relevant to diffractive J/ip photoproduction at HERA. 

Table 1: The correction factor F"^ of (23) evaluated using three recent sets of partons. 



X 


F2 (GRV) 


F^ (A') 


(G) 


5 X 10"^ 


0.46 


0.54 


0.48 


10-3 


0.45 


0.52 


0.47 


2 X 10-3 


0.45 


0.50 


0.47 


5 X 10-3 


0.44 


0.48 


0.46 



Besides the correction factor coming from the J/ip wave function weighted by the 7 — cc 
coupling, we have to consider the effect arising from the inequality 2mc 7^ M^. At first sight 
such a kinematic effect might appear to be negligible. However we note that eq. (2) is written 
in terms of M^, while in eqs. (21) and (22) the current quark mass rric enters. Thus we see that 
eq. (2) comes from eqs. (21) and (22), in the small limit, with replaced by |M|. The 
other factors, such as the leptonic width and the photon propagator, depend directly on M^. 
So for photoproduction, with = 0, the total relativistic correction factor is estimated to be 

.2 f U 



with a large uncertainty of at least ±30%. Thus it turns out that the best prescription is 
as follows. We should replace ml by in eqs. (21) and (22) so that the mass factor 

disappears from eq. (25) and the correction factor becomes closer to 1. Of course this 
does not decrease the overall uncertainty. However we emphasize that the uncertainty mainly 
affects the normalization of the gluon density extracted from J/ %p diffractive data, and that the 
prediction for the x dependence should be much more reliable. We see from Table 1 that the 
gluon density leads to a small x dependence of F^, which we neglect in our comparison with 
the data below, but which should be included when the J / ip data become more precise. 

We see from (21) and (^2]) that the suppression F^, arising from the use of the relativistic 



wave function, decreases with increasing Q due to the presence of the z{l — z)Q term. In fact 
since 

z{l-z) ^ \ {l-{vl)), (26) 

where {vD = we even expect, as Q'^ increases, that eventually the factor will give an 

enhancement, not a suppression, of the cross section. When we take Fermi motion into account 
the ratio of diffractive J/ip production from longitudinal and transverse photons becomes 



ctl Q 



2 



ax 4m^ 



(1 - {v')Y. (27) 
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Our result disagrees with the conclusions of ref.|ll] where a large suppression (of at least 
a factor of 3) was found to arise from Fermi motion. There are several differences in our cal- 
culations. First ref. |TT| uses a potential model J/ip wave function whereas here we employ a 
Gaussian form. Our considered range of (f ^) embraces a large class of wave functions and covers 
the results obtained using the wave function of any reasonable charmonium model. (Inciden- 
tally it is worth mentioning that the recent lattice calculation [|14| gives a wave function better 
described by a Gaussian form than by the more sophisticated potential models). However, 
most of the difference between our results and that of ref. [|l^ has a different origin. Part of the 
difference is due to the inclusion of the mass factor in (25). The remaining discrepancy arises 
from the use of the appropriate wave function of the virtual photon and the use of the correct 
scale for the gluon in eq. (21). 



(d) Rescattering or absorption of cc quark-pair 



The shadowing effects of the gluons are already included in xg{x, Q ) if the distribution 
is determined in a global parton analysis of experimental data. Here we are concerned with 
the rescattering or absorption of the cc pair as it transverses through the proton. A typical 
diagram is shown in Fig. 5. Let us study the first shadowing correction (with n = 1) arising from 
the exchange of one extra pair of gluons with transverse momenta -^Kt, where the Sudakov 
decomposition (see (|D) of the gluon 4-momentum is 















X 




" X 











Figure 5: A typical diagram leading to a cc rescattering correction. In the text we compute 
the n = \ contribution. All quark lines are on-shell except for that with momentum r, as in 
Fig. 3. 



= K = aKV \ pKq' + Kt- (28) 

To compute the extra gluon loop we first close the contour of the integrations over uk and (3k 
on the quark propagators marked by small crosses on Fig. 5. Hence we again have all quarks 
on- mass-shell, except for that with momentum r. Thus we do not change the result of the 
integration over the cc loop. If we take into account the four different ways the extra gluons 
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can couple to the c and c, then we have to include a total of 16 diagrams, and the integration 
over and gives, in analogy with (0), an amplitude of the form 



r (fkr 




1 ■ 


i krp J 




2 . 



<7^ 



^2 



^2 



+ 



1 



^2 



+ 



Q Q Q +\kT + KT 



/BFKL(a;, /ct) /bfkl(5', ^I) 



Q +\kT-KT\'' 



(29) 



We again neglect kt- The non-zero value of x arises because we have to put a quark on- 
shell with an additional momentum Kt- We therefore need to transfer to the quark a proton 
momentum fraction 



X 



X 



2zW^ 



(30) 



2 

since typically we have ~ Q /2. For kx, Kt <^ Q, after neglecting the angular correlation 
between kx and Kt, that is setting 



{{kT-KT)') = ^k^K^T, 



the expression in square brackets in (E9[) becomes 



1 
2 



4 4if| 



2 

We expand the integrand in (|29|) assuming that i^l^ -C Q . If we keep the leading terms 
then we find 



2 



3b 



dK^ 



^2 



d{xg{x,K^)) 
dK^ 



^T --'T 

where Q is defined so that the correction to the amplitude is 



A = Ao{i - in). 



(31) 



(32) 



As the correction is not too large we can eikonalize ( [31| ) by hand, directly in the momentum 
representation, and write the final expression in the form 



cr = ctq exp(— i7). 



(33) 



The factor 1/6 in (pT]) comes from the integral over the "reggeon" loop in Fig. 5 



dt' e 



-bt' 



(34) 



where we use the experimental slope 6 = 4.5 GeV ^. The Kt integration in (^) is evaluated 
just as in ([T7|) and (pTSl). 

It is worthwhile mentioning that the simple estimate of (31) is in good agreement with the 
more detailed treatment of ref. [0. 
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(e) Next-to-leading order corrections 



There are expectations that higher-order radiative corrections will have a negligible effect. 
We refer to four possible types of radiation shown all together in Fig. 6. 

The contribution of the diagram with the single extra gluon denoted by (a) is already 
incorporated in T{J/ip e"'"e"), while (b) is hidden in as and (c) is subsumed in xgi^x^Q"^). 
Only correction (d) survives, but this is suppressed by a factor Tint/T/, which is small for high 
W , that is for x ^ 1, see the comment in the caption to Fig. 1. 




Figure 6: Four different types of radiative correction, denoted by (a,b,c,d), are shown together 
on the lowest-order diagram for diffractive J/ ip production. 

Independent evidence for small radiative effects comes from a study of next-to-leading order 
corrections in inelastic J/ip production at HERA |T^. The lowest-order subprocess is shown in 
Fig. 7. It is found that the corrections to inelastic J/ip production are in general appreciable 
( > 50%). However, they reduce almost to zero when s M^, that is when the emitted gluon 
k is very soft, which corresponds to our pomeron with \x'\ <^ x, see (p!^). 

(f) Inclusion of the real part 

So far we have calculated only the imaginary part of the amplitude. At high energy W, 
that is at small x, our positive-signature amplitude behaves as 

A (X x~^ + {-xy^. 

Providing A is small, the amplitude is dominantly imaginary and the real part can be calculated 
as a perturbation 

ReA n n din A ti d\n{xg{x,Q )) . , 

Ml 2 2 dln{l/x) 2 dln{l/x) ' ^ ' 

12 



7 



s 




k 



Figure 7: The lowest-order subprocess, J I'^g-, describing inelastic J/V' production. The 

emitted gluon k is necessary to satisfy the colour and spin constraints of the J/V'. 

3. Comparison with diffractive J/?/; photoproduction data 

We use the perturbative QCD formula (|19D together with the corrections which we detailed 
in section 2, to calculate the cross section for diffractive J/'?/' photoproduction 



as a function of the 7p centre-of-mass energy W . We take the experimental value for the 
slope parameter: 6 = 4.5 GeV~^. The prediction depends on the square of the gluon density, 
\^g{x^Q )]^, at X = M'^jW'^ and for values of Q in the region of M^. For illustration we 
calculate the I^-dependence of the cross section using the gluon distributions of three of the 
latest sets of partons, namely GRV and MRS (A', G) The latter partons are extrapolated 
below = 4 GeV^ with next-to-leading GLAP evolution. 

It could be argued that it is inappropriate to use GLAP-based gluons in our resummed 
log(l/a;) formalism. However, they should be regarded simply as a parametrization of the 
data. The non-perturbative parameters are a good representation of all physical effects such 
as BFKL smearing and shadowing corrections. It is well known that the present data cannot 
distinguish the underlying perturbative physics. However, data do distinguish between the 
parametrizations. 

The three predictions are compared with the available high energy Jjt\} photoproduction 
data in Fig. 8. The curves correspond to the choice of lower limit Q\ = \ GeV^ in the integral 
over the gluon A;|. in ( |T9|) . We explored the sensitivity of the predictions to variation of Q\ over 
the range 0.5 to 2 GeV^ for the GRV gluon, and 1 to 4 GeV^ for the MRS gluons. Depending 
on the parton set, the cross section values change by ±15% at = 10 GeV, increasing to 
±25% at = 100 GeV as Ql is varied. 



1 da 
h Ht 
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Figure 8: The measurements [5,6] of the cross section for diffractive J/ip photoproduction 
compared with the full perturbative QCD prediction, as described in section 2, obtained from 
the three latest sets of partons. 

In Table 2 we quantify the effect of sequentially improving the prediction of the lowest- 
order formula (H), first to (|l^), which goes beyond the leading InQ^ approximation to include 
the effects of the transverse momenta of the exchanged gluons, and then to include the factor 
exp(— n), with Q given by (0), to allow for the cc rescattering on the proton. Table 2 shows 
the size of the effects at W = 100 GeV. We see that the two corrections to the lowest-order 
formula partially compensate each other. The corrections decrease with decreasing energy W, 
and are found to be insignificant for W ^20 GeV. 



Table 2: The cross section for diffractive J/ip photoproduction 



Partons 


(t(7P ^ J/ijp) 


in nb atW = 100 GeV 


lowest-order 


+ gluons kx 


-|- cc rescatt. 


GRV 


128 




157 


88 


MRS (A') 


47 




79 


55 


MRS (G) 


262 




331 


208 
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In addition to the above effects, the predictions shown in Fig. 8 also include the real part 
contribution via p5|). According to ( pSj) the J/ip relativistic effects are estimated to leave the 
cross section unaltered, but to introduce a sizeable normalization uncertainty. 

4. Conclusions 




Figure 9: The values of F2 versus In Q"^ aXx = 3.5 x 10"^ obtained from GRV [7] and MRS(A',G) 
[8] partons. The shaded band delimits the most recent measurements of the structure function 
obtained by the HI and ZEUS collaborations [18]. 

We conclude that the HERA J/ip photoproduction data offer a unique opportunity to 
distinguish between the various gluon distributions in a kinematic region (x = M'^/W'^ 10^^ 
and Q ~ 2.5 GeV^) where they are particularly distinct. In order to make a meaningful 
comparison we have computed several corrections to the lowest-order formula for the J/ip 
cross section. The uncertainties are found to be greater in the predicted size of the cross 
section than in the W or, equivalently, the x dependence. That is the shape, rather than the 
normalisation, is the better discriminator between the various gluons. The comparison of the 
J/ ip cross section data with the values calculated using the gluons of three of the latest parton 
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sets was shown in Fig. 8. The shape (and the normahsation) predicted by the MRS(A') gluon 
is, within the theoretical accuracy, in excellent agreement with the measured values of the J/'^ 
photoproduction cross section. In fact the data appear to favour the MRS (A') gluon over the 
GRV gluon, and to rule out that of MRS(G). 

A similar preference for the MRS (A') gluon is found by inspecting the slope of the most 
recent HERA measurements of the structure function F2 as a function of InQ^ as shown, for 
example, in Fig. 9. This slope is an independent measure of the gluon at small though, as 
anticipated, it is not such a sensitive probe as the W dependence of the HERA diffractive J /ip 
cross section data. However, if combined together, the HERA measurements of the diffractive 
cross section cr(7p — > J/ipp) and of the slope 9^2/(9 In offer an excellent opportunity to pin 
down the gluon density for x < 10~^. 
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